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La  Geometrie  de  Compensation  Non-Lineaire  - Le 
Probleme  Spatial  d’Intersection 
dans  l’Option  de  la  Geodesie  Tridimensionnelle 

Abdelmajid  Ben  Hadj  Salem,  Dipl.  Ing.0 


1 Introduction 

Dans  un  article  [Tj  E.  Grafarend  et  B.  Schaffrin  ont  etudie  la  geometrie  de 
la  compensation  ou  ajustement  non-lineaire  pour  le  cas  du  probleme  d’intersec- 
tion  plane  en  utilisant  le  modele  de  Gauss  Markov,  par  les  moindres  carres.  Le 
present  papier  developpe  la  meme  methode  en  travaillant  sur  un  exemple  de  la 
determination  d’un  point  par  trilateration  dans  l’option  de  la  geodesie  tridimen- 
sionnelle pour  la  determination  des  coordonnees  (x,  y,  z ) d’un  point  inconnu  a 
partir  des  mesures  des  distances  vers  n points  connus. 


2 La  Geometrie  Non  Lineaire  du  Modele  de  Gauss- 
Markov 

Le  modele  non  lineaire  de  Gauss-Markov  est  defini  par  : 

C(X)  = L-e;  e€AA(0,r)  (1) 

avec  : 

- L : le  vecteur  des  observations  (n  x 1)  = (Li,  L2, Ln)T , 

- X : le  vecteur  des  inconnues  (m  x 1)  = (X\,X2,  ..,Xm)T, 

- e : le  vecteur  des  erreurs  (nxl)  = (ei,  e2, en)T  suit  la  loi  normale  A^(0,  T) 
avec  E(e)  = 0 et  T = E(eeT)  la  matrice  de  dispersion  ou  variance,  on  prendra 
T = ctq.P”1.  P est  la  matrice  des  poids  et  do  une  constante  positive. 

- C : est  une  fonction  donnee  injective  d’un  ouvert  U C Mm  — >•  Mn  et  m < n. 

Remarque  : dans  le  cas  d’un  modele  lineaire,  la  fonction  ( = A.X  ou  A est 
une  matrice  n x m. 


On  note  Im(  = {(,{X)/X  £ U}  l’image  de  U par  la  fonction  £.  ImC,  est  une 
variete  de  dimension  m verifiant  les  conditions  : 


(i)  : les  vecteurs 
chaque  point  X € U, 


d(  d(  d( 

djT1,dx^,"',dx^ 


sont  lineairement  independants  en 


(ii)  : les  fonctions 


d2C 

dXidXj 


sont  continues  sur  U pour  i.j  £ {l,2,...,m}. 


1.  6,  rue  du  Nil,  Cite  Soliman  Er-Riadh,  8020  Soliman,  Tunisie 
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(2) 


On  introduit  un  produit  scalaire  : 

< Cl,  <2  >=  (T-P-C'2 

D’ou  la  norme  du  vecteur  £ = ((1,(2,  ■■■■,( n)T  ■ 

licn2=<cJo=cTJp-c  = E^-cf  (3) 

i= 1 

dans  l’espace  vectoriel  Mn  en  prenant  la  matrice  de  poids  P une  matrice  diago- 
nale  . 


Alors  la  solution  par  les  moindres  carres  X sera  definie  par  : 

\\L-((X)\\=min{\\L-aX)\\/X  eU}  (4) 

Cette  condition  est  exprimee  par  les  equations  de  Lagrange-Euler  soit  : 


d 

dx, 


\L  — ((X)\\2  = 0 pour  i G {1, 2, m} 


(5) 


En  effet,  on  veut  minimiser  la  fonction  : 

F(X)  = F(X1,X2,..,Xm)  = ||L-C(X)||  = \\L-C(X1,X2,...,Xm)\\  (6) 

Comrne  F est  une  fonction  positive,  minimiser  F c’est  aussi  minimiser  F2,  soit 
G(X)  = F2(X).  En  appliquant  les  equations  de  Lagrange-Euler,  on  obtient  : 


soit 


or  : 


d 

dXr 


dG(X) 
dX~ 

\L-C(X1,X2,...,X, 


= 0 


dG{X) 

dX, 


= 0 


= 0 pour  i € {1,  2, ...,  m} 


(7) 


I L - C(X1,X2, Xm)\\2  = {L-  C(XUX2,  ...,Xm)f.P.(L  - C(X1,X2,  ...,Xm))  = 

c (X)T.P.((X)  - 2Lt.P.((X)  + lt.p.l 

(8) 


Soit  : 


dG(X)  _n/_,v,T  DdC(X)  _2LT  p daX) 


dXi 
ou  encore  : 


= 2Q(Xy  .P.- 


dXi 


pour  i € {1,  2, ...,  m}  (9) 


dG(X)  t,t  „dC(X) 


dXi 


= 2 (aX)-Ly.P.- 


dXi 


dXi 

pour  i € {1, 2, m} 


ce  qui  donne  en  utilisant  © : 

d((X) 


<L-  C(X). 
ou  < e, 


dXi 

dgx) 

dXi 


>=  0 pour  i € {1, 2, m} 

(10) 

>=  0 pour  i € {1, 2, m} 

(11) 
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Geometriquement,  cela  veut  dire  que  le  vecteur  erreur  e = L — ( '(X ) est  per- 
pendiculaire  (produit  scalaire  nul)  au  plan  tangent  de  la  variete  ImC,  au  point 
C(X)  (s’il  existe). 


Pour  le  cas  non-lineaire,  la  condition  (1111)  est  necessaire  mais  non  suffisante. 

d2G 

Pour  obtenir  le  minimum,  il  faut  que  la  matrice  1 , — , , — , ),i,  j € {1,2, 


soit  definie  positive. 


3 Etude  d’un  cas  pratique 

On  considere  la  determination  d’un  point  par  trilateration  dans  l’option  de 
la  geodesie  tridimensionnelle  pour  la  determination  des  coordonnees  (x,  y,  z ) 
d’un  point  M inconnu  a partir  des  mesures  des  distances  vers  4 points  connus 
Mk(uk,  Vk,Wk)k= 1,2,3,4- 


Pour  faciliter  les  calculs,  on  prendra  P = I,  cro  = 1 et  nous  adoptons  la 
fonction  £(X)  comme  suit  : 


cw 


Ci  = (x  - Ui )2  + {y-  Vi)2  + (z-  Wx)2 

C2  = (x  - u2 )2  + (2/  - v2)2  + (z-  w2)2 

(3  = (x-  U3)2  + (y  - V3)2  + {z-  W3)2 

C4  = (x  - U4)2  + (y  - V4)2  + (z-  W4)2 


avec  X\  = x,  X2  = y,  X3  = z les  inconnues 


(12) 

(13) 


et  d’apres  (ED,  la  fonction  ( n’est  pas  une  fonction  lineaire  des  variables  X*.  ( 
est  une  fonction  de  M3  — >•  M4  qui  s’ecrit  : 


C(X)  = (kek 


ou  ek  est  la  base  orthonormee  de  Mn.  Voyons  qu’elle  verifie  les  deux  conditions 
(i)  et  (ii)  cites  ci-dessus. 


Calculons 


dc_ 

dXt  ’ 


on  a alors  : 


2(x  — u 1) 

2{y  - xi) 

2(z  — wi) 

ac 

2(x  - u2)  d( 

2 (y  - v2)  d( 

2 (z  — w2) 

dXx  " 

2(x  - U3)  ’ dX2 

2 {y  - v3)  ’ dX3 

2 (z  - w3) 

2(x  — U4) 

2 (y  - vA) 

2 (z  — W4) 

(14) 


Pourque  les  3 vecteurs  soient  lineairement  independants,  il  faut  que  les  points 
M,  Mi,  Mj , Mk  ne  soient  pas  alignes.  Pour  la  condition  (ii),  on  a facilement  : 


2 

2 

d2C 

2 d2( 

2 d2C 

dX2 

2 ’ ox2 

2 ’ dX2 

2 

2 

(15) 
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et  pour  i / j,  on  a : 


d2( 

OXidXj 


0 

0 

0 

0 


(16) 


Done  les  quantities 


d2gx) 

dXidXj 


sont  continues  et  la  condition  (ii)  est  verifiee. 


3.1  Ecriture  des  Equations  de  Lagrange-Euler 

Pour  determiner  la  solution  par  les  moindres  carres  du  rnodele  non-lineaire, 
on  ecrit  les  conditions  ([TUP.  Le  vecteur  L = {L\,  L2,  L3,  L4)r  telque  chacun  des 
Li  represente  le  carre  de  la  distance  spatiale  mesuree.  On  a alors  en  utilisant 

m ■ 

i = 1,2,3 

soit  : 

( x - ui)(L!  - Ci)  + (x-  u2)(L2  - C2)  + {x-  u3)(L3  - C3)  + {x  - u4)(L4  - C4)  = 0 

(y  - ^i)(^i  - Ci)  + (y  - v2){L2  - c2)  + {y  - v3)(L3  - ( 3 ) + (y  - u4)(L4  - c4)  = 0 

(z  - wi)(Li  - Cl)  + (z-  w2)(L2  - C2)  + {z  - w3)(L3  - C3)  + (z-  w4)(L4  - C4)  = 0 

(17) 

Les  equations  (1171)  represente  un  systeme  de  trois  equations  non  lineaires  de 
trois  inconnues  (x,  y,  z)  dont  la  solution  est  un  peu  compliquee. 

Pour  faciliter  encore  la  resolution  du  systeme  precedent,  on  va  supposer  que 
la  variable  2 est  connue  egale  a zo,  dans  ce  cas,  on  se  limite  a trois  distances 
mesurees  L\,L2  et  L2.  Alors  (1171)  s’ecrit  : 

(x  - ui){L\  - Ci)  + (x-  u2)(L2  - C2)  + (x-  u3)(L3  - C3)  = 0 
(y  - vi)(Li  - Ci)  + (y-  v2 ){L2  - C2)  + (y-  v3){L3  - C3)  = 0 (18) 


avec  : 

Ci  = (x  - Ui)2  + (y  - U1)2  + (z0  - Wi)2 

C2  = {x  - u2)2  + (y-  v2)2  + (2:0  - w2)2 

C3  = (x  - u3)2  + (y-  v3)2  + ( z0  - w3)2 

Les  expressions  Ci  — L,  s’ecrivent  sous  la  forme  : 

(i  - Li  = x2  + y2  - 2 xui  - 2yvi  + a,  (19) 

avec  a*  = constante  (20) 
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Le  systeme  (1151)  devient  : 

(x  - u\)(x 2 + y2  - 2xui  - 2yv\  + ai)  + (x  — u2)(x2  + y2  - 2 xu2  - 2 yv2  + a2) 

+(x  - u3)(x2  + y2  - 2xu3  - 2yv3  + a3)  = 0 

(21) 

(y  - vi)(x2  + y2  - 2xu\  - 2yv\  + ai)  + (y  — v2)(x2  + y2  - 2xu2  - 2 yv2  + a2) 

+(y  - v3)(x2  + y2  - 2 xu3  - 2 yv3  + a3)  = 0 

(22) 


En  developpant  les  equations  Oil)  et(021).  on  obtient  : 

3x3  + 3 xy2  — 3x2(ui  + u2  + tt3)  — y2(iti  + u2  + u3)  - 2xy(ui  + u2  + u3)+ 

x(ai  + a2  + a3  + 2 u\  + 2u2  + 2 u\)  + 2y(u\V\  + u2v2  + ^3^3)  — (a\Ui  + a2u2  + a3u3 ) = 0 

(23) 

3 y3  + 3 yx2  — 3y2(tq  + v2  + f3)  — x2(ui  + u2  + u3)  - 2xy(u\  + u2  + u3)+ 

y(ai  + a2  + a3  + 2u2  + 2v\  + 2u2)  + 2x(u\Vi  + u2v2  + u3v3)  — (aitq  + a2v2  + a3v3 ) = 0 

(24) 


Supposons  qu’on  se  limite  a deux  distances  L\  et  L2,  alors  on  a a resoudre  : 

2x3  + 2 xy2  — 2x2(ui  + u2)  - y2{u\  + u2)  — 2xy(v\  + v2)+ 
x(ai  + a2  + 2 u\  + 2u2)  + 2y{u\V\  + u2v2)  - (aiui  + a2u2 ) = 0 (25) 

2 y3  + 2 yx2  — 2y2(v\  + v2)  — a:2(ui  + v2)  — 2xy(u\  + u2)+ 
y(ai  + a2  + 2u2  + 2u|)  + 2x(u\V\  + u2v2)  — (a^i  + a2v2)  = 0 (26) 


3.2  Reduction  des  Equations  de  Lagrange-Euler 

Dans  ce  paragraphe,  on  essaye  de  reduire  l’ecrirure  du  systeme  (fOUl)  - f 1 2 6 . 
A cet  effet  posons  : 


z = x + iy,  z = x — iy  avec  i = \/—l 
s = ui+u2,  t = v i + v2,  p = uf  + ul,  q = v\  + V2 
a = ai+a2,  r = u\V\  + u2v2,  d = a\U\  + a2u2,  f = a\vi  + a2v2  (27) 

ce  qui  donne  : 

2x  = z + z,  2 iy  = z — z (28) 

Alors  les  expressions  (12511  - (12611  deviennent  : 


4 z2z  + Azz2  + (2  it  — s)z 2 — 6 szz  — (s  + 2 it)z2  + 2 z(a  + 2 p — 2 ir)  + 2z(a  + 2 p + 2 ir)  — Ad  = 0 

(29) 

4 z2z  — 4 zz2  + (it  — 2 s)z2  — 6 itzz  + (2s  + it)z2  + 2 z(a  + 2 q + 2 ir)  — 2 z(a  + 2q  — 2 ir)  — 4if  = 0 

(30) 


Nous  presentons  dans  la  suite  la  resolution  des  equations  ci-dessus  : 
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4 Resolution  du  Systeme 


Posons  : 


jO 

z = x + iy  = pe 
s + lit  = lelu 
a + 2p  + 2 ir  = me” 
o + 2q  + 2ir  = ke 
2s  + it  = helip 


(31) 

(32) 

(33) 

(34) 

(35) 


Utilisant  les  relations  : 


z = 2 ipsind, 


z + z = 2 pcosO 
sind,  zz  = p2 


(36) 

(37) 


Les  equations  (T2U1)  et  (|5U1)  deviennent  : 

8 p3cos0  — 6 sp2  — 2lp2cos{29  — u)  + 4 mpcos(Q  — a)  — M = 0 (38) 

8 ip3  sind  — 6itp2  — 2ihp2  sin{26  — </?)  + 4 ikpsin(0  + p)  — Aif  = 0 (39) 

Soit  : 

4 p3cosd  — [3s  + lcos{26  — w)]  p2  + 2 mpcos(6  — a)  — 2d  = 0 (40) 

4 p3sind  — [3i  + hsin{26  — p)\  p2  + 2 kpsin(6  + /r)  — 2/  = 0 (41) 

Eliminant  les  terrnes  constants  entre  (I40p  et(!41l)  et  apres  simplification  par 
p yA  0,  on  obtient  l’equation  : 

4 p2(dsind  — fcosO)  + p(3sf  — 3 td  + lfcos(20  — oj)  — hdsin(20  — </?))+ 


2kdsin{6  + p)  — 2mfcos{6  — a)  = 0 (42) 


Maintenant,  on  elimine  les  coefficients  en  p3  des  equations  ([4UD  et()4ip.  nous 
obtenons  une  equation  en  deuxieme  degre  en  p : 

p2(3s.sin9  — 3tcos9  + lsin6cos(26  — oj)  — hcos6sin{29  — </?))+ 

2 p{kcos9sin{9  + p)  — msin9cos(9  — a))  — 2 fcos9  + 2 dsin9  = 0 (43) 

5 Resolution  du  systeme  des  inconnues  en  p et  9 

On  ecrit  les  equations  (EP1)  et(|43p  sous  la  forme  : 


Ap2  + Bp  + C = 0 
A'p 2 + B'p  + C'  = 0 


(44) 

(45) 


Le  systeme  precedent  est  resolvable  si  et  seulement  si  : 


A _ B _ c 
A'  ~ ~B'  “ ~C' 
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(46) 


soit  en  choisissant  : 


(47) 


§ = §=*•  BC’-B'C  = 0 


Soit  : 


[3s  f — 3 td  + lfcos(29  — oj)  — hdsin(29  — tp)\  (2 dsinO  — 2fcos6)  = 

[2 kcos9sin(9  + y)  — 2 msin9cos(9  — ct)]  [2fca!sin(0  + //)  — 2 mfcos(9  — a)] 

(48) 

ce  qui  donne  un  polynome  en  sind,  cosd,  sin29,  cos29.  En  posant  : 

£ = tg9  (49) 

et  utilisant  les  formules  : 

2£ 


sin29  = 
cos29  = 


1 + £2 
i-£2 
1+  £2 


(50) 

(51) 


L’equation  devient  un  polynome  du  troisieme  degre  en  ^ qu’on  peut  re- 
soudre  par  les  methodes  classiques. 


6 Calcul  de  la  Matrice  Covariance  des  inconnues 

Comme  on  a : 

Ci  = {x  - ui)2  + {y  - vi)2  + (z0  - wif 
C = I C2  = (x  - u2 )2  + (y-  v2 )2  + (zo  - W2)2 
C,3  = {x-  u3)2  + (y-  v3)2  + (z0  - w3) 2 


(52) 


Or,  on  s’est  limite  a deux  inconnues,  le  vecteur  £ = (Ci,(2)T-  D’ou  en  differen- 
tiant  (1521)  en  utilisant  les  deux  premieres  lignes,  on  obtient  alors  : 


d(  = 


d Ci 

d(  2 


2(x-ui)  2{y-vi) 

2(x  - u2)  2 (y  - v2) 


dx 

dy 


avec  : 


Jx  = 2. 


x — u\  y — v 1 
x — u2  y — v2 


= Jx-dX  (53) 


(54) 


Par  suite 


dX  = 


dx 

dy 


= J 


-1 
x ■ 


dCi 
d(  2 


d’ou  la  matrice  covariance  du  vecteur  inconnu  X = (x,  y)1 

2 T -1  2 T -iT 
&x  = dx  -V L-dx 


(55) 


(56) 
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